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Abstract. The NA60 experiment has studied low-mass muon pairs production in proton- 
nucleus collisions at 400 GeV/c at the CERN SPS. The mass spectrum is well described by the 
superposition of the two-body and Dalitz decays of the light neutral mesons 77, p, cj, 7/ and <j). 
The large collected statistics allows to isolate the contributions due to the Dalitz decays of the 
77 and uj mesons, from which the electromagnetic transition form factors of the two mesons have 
been extracted. The found values agree with the previous available measurements, improving 
their uncertainty thanks to the higher statistics. The results thus confirm the discrepancy with 
the prediction of the Vector Meson Dominance model in the case of the electromagnetic form 
factor of the uj meson. 



Proton-nucleus studies have been part of the heavy-ion experimental program at the CERN 
SPS, of which the NA60 experiment represents the most recent and one of the most fruitful 
efforts. Beyond serving as a robust reference baseline for the heavy-ion data, proton-nucleus 
data also stimulate a genuine physics interest, allowing measurements in an environment of 
cold nuclear matter. This is in fact the case of the electromagnetic transition form factors of 
the 77 and uj mesons, which have been studied here through their Dalitz decays 77 /i + /i~7 
and uj —± /i + /i - 7r°. In these A — >► Bl + l~ processes, the mesons decay electromagnetically into 
a virtual photon with mass M — in turn converting into a lepton pair — and a third body. 
Assuming point-like particles, the mass distribution of the emitted lepton pair (that is, the 
decay rate of the process as a function of M) can be exactly described by QED as firstly shown 
by N.M. Kroll and W. Wada [I] in the 1950s. 

However, the observed mass distribution for the lepton pairs emitted in an A — >> Bl + l~ decay 
differs from the point-like one, because of the internal electromagnetic structure of the decaying 
meson, which reveals at the vertex of the transition A — >> B. This difference can be formally 
described by a function of M, usually expressed as \Fab(M)\ 2 , known as transition form factor. 
Experimentally, \Fab(M)\ 2 is directly accessible by comparing the measured invariant mass 
spectrum of the lepton pairs from the Dalitz decays, with the point-like QED prediction. The 
transition form factors are usually given in terms of the pole parametrization: 

|F| 2 = (1-M 2 /A 2 )- 2 . (1) 

The Vector Meson Dominance (VMD) model gives theoretical predictions both for A~ 2 and 
A~ 2 . These predictions have been tested by the Lepton-G experiment [3] using pion beams, and 



recently by the NA60 experiment in In-In peripheral collisions [3]. While the experimental values 
agree among each other both for the 77 and uj mesons, only the VMD model prediction for the 77 
meson is confirmed by the data, while the one for the uj meson is not, strongly underestimating 
the observed form factor. The NA60 p-A data can now give a further significant contribution 
to the measurement, thanks to the high quality and the large statistics available. 

A description of the NA60 apparatus can be found for example in [5j. Here, only few relevant 
details are given. The produced dimuons are identified and measured by the muon spectrometer, 
composed of a set of tracking stations, trigger scintillator hodoscopes, a toroidal magnet and 
a hadron absorber. The stopping of hadrons (provided by more than 20 interactions lengths) 
allows to select the highly rare dimuon events. On the other hand, the material which stops 
the hadrons also induces multiple scattering and energy loss on the muons, degrading the mass 
resolution of the measurement made in the spectrometer. To overcome this problem, NA60 
already measures the muons before the absorber with a tracking telescope, made of pixel silicon 
detectors jB]. The muon tracks are found among the many other charged particle tracks by a 
correct matching, both in angular and momentum coordinates, with the reconstructed tracks in 
the muon spectrometer. 

The preliminary results presented here are obtained from the analysis of the data collected 
during three days within the 2004 proton run, with a system of nine sub-targets of different 
nuclear species, simultaneously exposed to an incident 400 GeV proton beam. The first step 
of the analysis is a fit of the reconstructed mass spectrum (see Figure [I]). The fit is performed 
comparing the data with the superposition of the expected sources (estimated via the MC). 
These sources are the Dalitz decays of the 77, uj and 7/ mesons, the 2-body decays of the 77, p, 
uj and (f) and the contribution of the open charm process. Before the comparison between data 
and MC is performed, the small component of combinatorial background (originating from tt 
and K decays) is subtracted from the real data: its shape is estimated with an event mixing 
technique, while its normalization is established fixing the LS (like-sign) component coming 
out from the mixing, to the LS component of the data (containing no signal from correlated 
pairs at the SPS energies). The background accounts for less than 10% of the integrated mass 
spectrum below 1.4 GeV/c 2 . The comparison between mixed and real sample gives an average 

uncertainty of 10%, for both the (++) and the ( ) components: given the rather low level of 

the combinatorial background contribution, this uncertainty hardly affects the evaluation of the 
shape for the Dalitz decays under study. 

The processes which are not of interest for the measurement of the form factors of the 77 
and uj are then subtracted from the real data, using the normalizations fixed by the fit. Only 
the Dalitz decays of 77 and uj are retained, together with the two-body decay of the p: this 
latter is left in, making it much easier to control the systematics due to the small contribution 
coming from the low-mass tail of the p in the mass region of interest here. The mass spectrum 
resulting after the subtraction is then corrected for the effects of geometrical acceptance and 
reconstruction efficiency. In order to do this, a correction profile as a function of the mass is 
built, weighting the profiles obtained from the MC for the three processes separately. After the 
correction for the acceptance x efficiency effects, the resulting mass spectrum — now free from 
any acceptance distortion — is compared with the superposition of the expected line shapes of 
the three processes 77 —> /i + /i~7 , uj —> p + p~7T° and p — >> p + p~ . 

Figure [2] shows the fit on the acceptance-corrected mass spectrum (black triangles), having 
isolated the contribution of the three processes 77 — >> p + p~j , uj p + p~7T° and p p + p~ . 
The fit function is a superposition of the expected line shapes for the three processes, described 
by the following functions: 
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In the fit, the following parameters are left free: 

• the three normalizations, one for each of the line shapes above 

• the parameters A~ 2 and Aj 2 , contained in the form factors \F T] (M 2 )\ 2 and |F W (M 2 )| 2 . 

Concerning the line shape for the p — » p + p~ decay, the parameters T p (entering in the definition 
of r to t) 5 Tn p and the effective temperature T p are fixed respectively to T p — 151 MeV (from 
the PDG), m p = 768 MeV/c 2 (also from the PDG) and T p = 170 MeV (consistent with the 
expected value obtained by statistical model fits on particle ratios in pp interactions). All 
the parameters involved in the analysis procedure have been varied to study their influence 
on the fits and determine the systematic uncertainties on the results: the weights defining the 
acceptance x efficiency correction, the fit range (and the choice of the free parameters), the value 
for the (Jrf/au ratio, the charge correlation parameter for the combinatorial background and 
finally the level of the open charm contribution. The analysis leads to the following values: 

A" 2 = 1.946 ±0.096 (stat.) ±0.015 (syst.) (GeV/c 2 )" 2 
A~ 2 = 2.248 ±0.030 (stat.) ±0.009 (syst.) (GeV/c 2 )" 2 

Once the final fit parameters and their errors are fixed, the contributions of the 77 p + p~~f 
and cj — )> p + p~7v° processes are disentangled, making it then possible to present the two form 
factors in the usual way, in order to compare them with the previous experimental points, see 
Figure |3j The pole parameters and their errors as obtained from the combined fit to both Dalitz 
decays are shown as inserts. Both figures include the NA60 data obtained in peripheral In-In, 
the Lepton-G data, and the expectations from VMD for comparison. Within the errors, perfect 
agreement between the three data sets is seen in both cases. Irrespective of the much reduced 
errors, the form factor of the 77 is still close to the expectation from VMD. The form factor of 
the uj 1 on the other hand, strongly deviates from the VMD, showing a relative increase close to 
the kinematic cut-off by a factor of ~ 10. 

In conclusion, we have presented a new measurement of the electromagnetic transition form 
factors of the 77 and uj mesons. The analysis has been performed by studying the Dalitz decays 
77 — >> /i + /i~7 and uo —> p + p~7T° in p-A collisions at 400 GeV/c with the NA60 apparatus. The new 
preliminary values presented here are in good agreement with those already available; improving 
the precision of the measure by a factor ~ 2, they confirm on more solid grounds the discrepancy 
between the VMD prediction for the uj and the experimental observations. 
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Figure 1. Fit on the target-integrated mass 
spectrum. 
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Figure 2. Fit on the acceptance-corrected 
mass spectrum. The shaded areas indicate 
the Kroll-Wada (KW) expectations for 
the r] and uj Dalitz decays for point-like 
particles, defined by QED [lj. 
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Figure 3. Electromagnetic transition form factors for the 77 (left) and uj (right) mesons as a 
function of M. 
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